Realistic modelling of quantum point contacts subject to high magnetic fields and 
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The electron and current density distributions in the close proximity of quantum point contacts 
(QPCs) are investigated. A three dimensional Poisson equation is solved self-consistently to obtain 
the electron density and potential profile in the absence of an external magnetic field for gate and 
etching defined devices. We observe the surface charges and their apparent effect on the confinement 
potential, when considering the (deeply) etched QPCs. In the presence of an external magnetic 
field, we investigate the formation of the incompressible strips and their influence on the current 
distribution both in the linear response and out of linear response regime. A spatial asymmetry of 
the current carrying incompressible strips, induced by the large source drain voltages, is reported 
for such devices in the non-linear regime. 
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INTRODUCTION 



The new era of quantum information processing at- 
tracts an increasing interest to investigate the intrinsic 
properties of small-scale electronic devices. One of the 
most interesting of such devices is the so called quan- 
tum point contacts (QPCs), where a quantized cur rent 
is transmitted through it under certain conditions 1 ^. 
They are constructed on two dimensional electron sys- 
tems (2DES) either by inducing electrostatic potential 
on the plane of 2DES and/or by chemically etching the 
structure. The essential physics is that, the small size 
of the constraint creates quantized energy levels in one 
dimension (perpendicular to the current direction) there- 
fore, transport takes place depending on whether the en- 
ergy of the electron coincides with this quantized energy 
or not. In the ideal case at low bias voltages, if the energy 
of the electron is smaller then the lowest cigcn-cncrgy of 
the constraint, no current can pass through the QPC. 
Otherwise, only a certain integer number of levels (chan- 
nels) are involved, therefore conductance is quantized 3 . 
Beyond being a useful play ground for the basic quantum 
mechanical applications many other interesting features 
are repo rted in the literature such as the 0.7 conductance 
anomal y! 4 ! 5 ^ 6 !, which became a paradigm since then. An- 
other adjustable parameter which induces quantization 
on the 2DES is the magnetic field B applied perpen- 
dicularly to the system. Such an external field changes 
not only the density of states (DOS) profile of the 2DES 
but also the screening properties of the system drasti- 
cally. The interesting physics dictated by this quanti- 
zation is observed as the quantum Hall effect iP™. Re- 
cent theoretical investigations point out the importance 
of the electron-electron interactions in explaining the in- 
teger quantum Hall effectP 1 ^, believed to be irrelevant 
in the early days of this fieloP"H which we discuss briefly 
in this work. The basic idea behind the inclusion of the 
interaction is as follows: due to the perpendicular mag- 



netic field the energy spectrum is discrete, known as the 
Landau levels (LLs), and is given by E n = hio c {n + 1/2) 
where n is a positive integer and cyclotron energy is de- 
fined as huj c = heB/m*, with effective electron mass 
m* = 0.067 * m e . Taking into account the finite size 
of the sample, i.e. the confinement potential, and the 
mutual interaction (Hartree) potential, one obtains the 
total potential. In the next step for a fixed average elec- 
tron density one calculates the resulting electron density 
distribution and from this distribution re-calculates the 
potential distribution iteratively. This self-consistent cal- 
culation ends in formation of the compressible and incom- 
pressible regions. In a situation where Fermi energy Ep 
is pinned to one of the LLs, then the system is compress- 
ible. Otherwise, if Ep falls in between two consecutive 
LLs, the system is known to be incompressible and, since 
there are no available states at the Ep for electrons to be 
redistributed, screening is poor. However, within these 
incompressible regions the resistivity vanishes, hence all 
the applied current is confined to these regions. We will 



be discussing the details of this model in Sec IV 



Most recently, the experiments performed at the 2DES 
including the QPCs, in the presence of an exter nal mag- 
netic field, manifested peculiar results^* 13 * 14 ! . In the 
first set of experi ments electron interference, such as 
Mach-Zehnder (MZ) 12 ' 15 ' 16 ' 17 ' 18 ' 19 ' and Aharonov-Bohm 
(AB^PHl, was investigated. The MZ interference experi- 
ments exhibit a novel and yet unexplained behavior, re- 
garding the interference contrast (visibility) at the inter- 
ferometer in the nonlinear transport regime (finite trans- 
port voltage). As a function of voltage, the visibility 
displayed oscillations whose period was found to be in- 
dependent of the path lengths of the interferometer, in 
striking contrast to any straightforward theoretical model 
(e.g. using Landauer-Biittiker edge states (LBES^PS). In 
particular, a new energy scale, of order of /ieV, emerged, 
determining the periodicity of the pattern. This unex- 
pected behavior of interfering electrons is believed to be 
related to electron-electron interactions^. The most sat- 
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isfying model existing in the literature is proposed by I. 
P. Levkivskyi and E. V. Sukhoruko-v 22 . However, other 
schemes also including interactions are presenP^ and 
models, which consider a non-Gaussian noise as a source 
of the visibility oscillations^, without interaction. The 
novel magnetic focusing experiments concerning QPCs 
has revealed the scattering processes taking place near 
these devices^. It was explicitly shown that, the exper- 
imental realization of the sample and the device itself 
strongly effects the transport properties. It is reported 
that, the potential profile generated by the donors (im- 
purities) and the gates deviates strongly from the ideal 
"point" contact. Even in each cool down process, since 
the impurity distribution changes, the quantum interfer- 
ence fringes differ considerably. Hence a realistic mod- 
elling of a QPC is desirable, which we partially attack in 
this work. Another interesting set of experiments within 
the integer quantum Hall regime is conducted by S. Rod- 
daro et.aP^, where the transmission is investigated as a 
function of the gate bias. The findings show that current 
transmission strongly deviates from the expected chiral 
Luttinger liquid^ behavior, since the transmission is ei- 
ther enhanced or suppressed by changing the gate bias. 
This effect was attributed to the particle hole symme- 
try of the Luttinger liquid and is discussed in detail in 
Ref. [25]. However, the explicit treatment of the QPC 
was left unresolved. Since the essential physics can be 
still governed by considering a finite size QPC opening, 
therefore assuming formation of a (integer) filling region 
sufficient. 

The theoretical investigation of QPCs covers a wide va- 
riety of approaches, which can be grouped into two: i) the 
models that include electron-electron interactions and ii) 
the models that do not. At the very simple model in de- 
scribing the QPCs, one considers a potential barrier per- 
pendicular to the current direction quantizing the energy 
levels. Therefore the electrons are considered to be plane 
waves before they reach to the QPC and transmission 
and reflection coefficients are calculated from this poten- 
tial profile. A better (2D) approximation is to model the 
QPCs with well defined function d 27 * 28 !, such as ellipses 
which lead to analytic solutions for the energy eigenfunc- 
tions and energies. About a decade ago J. Davies and 
his co-workers developed the "frozen charge" model 29 
to calculate the potential profile induced by the gates 
defining the QPC. This approach is still one of the most 
used technique to obtain the potential profile, however, 
it is not self-consistent and completely ignores the donors 
and surface charges. There exists many theories which 
takes the potential profile from the frozen charge model 
as an initial condition, and provides explicit calculation 
scheme s to obtain charge, current and potential distri- 
butiong^°- | 31 | 32 | 33 | 34 l One of the most complete scheme, 
even in the presence of an external B field, is the local 
spin density approximation (LSDA) within the den sity 
functional theory (DFT^P. The LSDA+DF T 31 ' 36 ' 37 ' ap- 
proaches are powerful to describe the essential physics of 
density distribution and even 0.7 anomaly phenomeno- 



logically, however, the description of the current dis- 
tribution is still under debate. The scattering prob- 
lem through the QPCs is usually handled by the "wave 
packet" formalism and is very successful in explaining 
the magnetic focusing experiments. However, the poten- 
tial profile is not calculated sclf-consistently and there- 
fore, the effect of the incompressible strips resulting from 
electron-electron interaction is not taken in to account. 

Back to early days of the theories that account for elec- 
tron interactions, i.e. Chklovskii, Shklovskii and Glaz- 
marPSl (CSG) and Chklovskii, Matveev and Shklovskii 27 
(CMS) models, the influence of the formation of the in- 
compressible strips has been highlighted. In the CMS 
paper, it was even conjectured that, 'the ballistic con- 
ductance of the QPC in strong magnetic field is given by 
the filling factor at the saddle point of the electron den- 
sity distribution multiplied by e 2 /2irh\ which is quantized 
only if an incompressible strip (region) resides at the sad- 
dle point. In one of the recent approaches, in the pres- 
ence of a strong B field, the electron-electron interaction 
is treated explicitly within the Thomas-Fermi approx- 
imation (TEA) self-consistently, meanwhile the current 
distribution is left unresolved^. In this model, similar 
to other approaches, the bare confinement potential is 
obtained from the "frozen charge" approximation, which 
in turn lead to discrepancies due to its non-self-consistent 
approach. Here, we improve on this previous work in two 
main aspects: i) the electrostatic potential is obtained 
self-consistently in 3D, which allows us to treat also the 
etched structures ii) the current distribution is calculated 
explicitly using the local version of the Ohm's law, also 
in the out of linear response regime. We organize our 
work as follows: In Sec|IT] we briefly describe the numer- 
ical scheme to calculate the potential profile at B = 
following Ref.JSS], which is based on iterative solution of 
the Poisson equation in 3D. In particular, we study the 
effect of different gate geometries and focus on the com- 
parison of the potential profiles of gate and etch defined 
QPCs. The numerical scheme to calculate potential and 
density profiles at finite temperature and magnetic field 
is introduced in Sec |III| Here, we review the essential in- 
gredients of the TFA and discuss the limitations of our 



approach. Sec IV is dedicated in investi gating t he current 
distribution within the local Ohm's law 9 * 40 * \ where we 
consider both the linear response (LR) and out of linear 
response (OLR) regimes. In the OLR, we show that the 
large current bias induces an asymmetric distribution of 
the incompressible strips, due to the tilting of the Lan- 
dau levels resulting from the position dependent chemical 
potential. We conclude our work by a summary. 



II. 



ELECTROSTATICS IN 3D 



The realistic modelling of 2DES relies on solving the 
3D Poisson equation for given boundary conditions, set 
by the hetero-structure (GaAs/AlGaAs in our calcula- 
tions) and the gate pattern, which describes the charge 



3 



and potential distribution. The hetero- structure, shown 
in Fig(l^,, consists of (metallic) surface gates (dark semi- 
elliptic regions on surface), a thin donor layer (denoted 
by light thin layer and S Silicon doping) which provides 
electrons to the 2DES and the 2DES itself indicated by 
minus signs confined to a thin area. The 2DES is formed 
at the interface of the hetero-structure. The average elec- 
tron density Wei (and its spatial distribution n e i(x,y)) is 
dictated by the donor density no and the metallic gates. 
Once the number of donors and the gate voltage Vg are 
fixed, the potential and charge distribution of the sys- 
tem can be obtained by solving the Poisson equation, 
self-consistcntly. 

For typical nanoscale devices with many (or at least 
a few) electrons in each of the electrostatically-defined 
regions, the charge distribution and the major energy 
scales are described to a good approximation by classical 
electrostatics. Due to the strong electric fields generated 
by segregating charge in a 2DES, the Coulomb energy 
is the dominant energy scale. In this sense, it is desir- 
able to have a self-consistent electrostatic description of 
the system if one expects a good quantitative description 
thereof. 

For solving the electrostatics of the system in three di- 
mensions we used a cod e dev eloped and successfully ap- 
plied in previous studies^Ml. it is based on a 4 th order 
algorithm operating on a square grid. The code allows 
flexible implementation of many boundary conditions rel- 
evant for nanoscale electrostatics: standard boundaries 
such as conducting regions at constant voltage (poten- 
tial gates), of constant charge (large quantum dots) or 
charge density (doping), but also boundaries such as a 
depletable 2DEG, dielectric boundaries and surfaces of 
semiconductors with the Fermi energy pinned due to sur- 
face charges. Since the calculation is constrained to a 
finite volume of space including the surface of the sam- 
ple, the outer boundary is considered open and is also 
obtained self-consistently along with the rest of the cal- 
culation. 

Overall the code provides a reliable description of the 
potential landscape and thus the electric field as well as 
the charge distribution for the sample under considera- 
tion. 

As an illustrating example in Fig(l]we show the hetero- 
structure under investigation together with the charge 
distributions at different layers. Area of the unit cell is 
1.5 x 1.5 /im 2 , whereas the hight is chosen to be 156 
ran. The donor and the electron layer lies 43 nm and 56 
nm below the surface, respectively. The metallic gates 
are deposited on the surface of the structure and are 
biased with —1.7 V and a homogeneous donor distribu- 
tion is assumed, Fig.[l];. The induced charge distribution 
on the metallic gates exhibits apparent inhomogencitics 
shown in Fig[TJj. The electrons are accumulated near the 
boundaries of the gates, a well known behavior for metal- 
lic boundary conditions. However, here we provide the 
explicit distribution which strongly differs from the pre- 
viously used " frozen charge" model where only a constant 
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FIG. 1: (Color online) (a) The Silicon doped hetero- 
structure, the 2DES is formed at the interface of the 
GaAs/AlGaAs (denoted by minus signs) and the metallic 
gates are deposited on the surface. At zero gate bias, the elec- 
tron density is determined by the number of donors, which we 
chose to be 4 x 10 16 cm -2 . Charge distribution at different 
layers, at the gates (b), the dopant layer (c) and the 2DES (d). 
It is clearly seen that not all the excess electrons are captured 
by the 2DES, rather a significant amount is accumulated on 
the surface. The electrostatic quantities are normalized with 
the relevant scales, i.e. charge (density) is normalized with 
the average electron density (e.g. Q2des(z, y) = n e i(x, y)/n cl ) 
and electrostatic potential (energy) with the potential energy 
of a single electron. 



potential profile is assumed. The influence of these in- 
duced charges become more important when considering 
an external B field. Since, the steepness of the external 
potential profile determines the effective widths of the 
current carrying incompressible strips. We should note 
that, our self-consistent model enables us also to handle 
the (side) surface charges which becomes important when 
considering chemical etching. In the following part we 
investigate systematically, the effects of the gate voltage 
and the device geometry on the electrostatic quantities. 
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FIG. 2: (Color online) Spatial distribution of the electrons 
as a function of the gate voltage. At zero bias (a) more elec- 
trons are populated under the gate which changes till deple- 
tion starts (b-d), where the gate voltage is set to be (b) -0.3V 
(c) -0.7V (d) -1.0 V. Almost no electrons are left beyond -1.5V 
applied to the gates, (e) -1.5V (f) -2.0 V. 



A. Gate defined QPCs 

In this subsection we compare the electron density pro- 
files calculated for different QPC geometries applying 
various bias voltages, which exhibits strong non-linear 
behavior in contrast to many models used in the litera- 
ture. We start our discussion with a rather smooth con- 
figuration, where the distance between the gates (W) is 
chosen to be 200 nm (see Fig{l]i). In Fig. [2] we show 
the cross section of the electron density profile for dif- 
ferent gate voltages. Interestingly, at Vg — we see 
that more electrons are residing beneath the gates. This 
effect is due to inhomogeneous (induced) charge distri- 
bution at the metallic gates similar to the distribution 
shown in FigjTJi. The induced charges are mostly ac- 
cumulated near the gate boundaries, whereas the center 
of the gates has almost a constant charge profile. In- 
creasing Vq to -0.3 Volts, already starts to depopulate 
electrons under the gates and the depopulation rate re- 



mains linear to the applied gate potential until the 2DES 
becomes depleted. In the [—0.5,-1.2] Volt interval, the 
density distribution changes relatively smooth, since the 
electrons can still screen the external potential quite well. 
It is important to recall that, in the absence of an ex- 
ternal B field, the DOS of a 2DES is just a constant 
D a (= m*/nh 2 = 2.83 x 10 10 meV" 1 cnr 2 for GaAs), 
which is set by the sample properties, therefore screen- 
ing is nearly perfect. Whereas this changes considerably 
when the electrons are depleted under the gate. This is 
observed by the strong drop of the potential when the 
depletion bias is reached around Vg — — 1.3 V. A sudden 
variation appears at the potential profile when the elec- 
trons are depleted beneath the gates (y — 350 nm) where 
the Vg becomes larger than -1.5 V. Therefore, the simple 
picture describing the QPCs as a smooth function of the 
applied gate voltage fails. In that picture it is assumed 
that the Fermi energy of the system remains constant and 
the potential profile, given by a well defined function, of 
the constriction is simply shifted by the amount of po- 
tential applied to the gates. Such a model is reasonable 
in the regime where the gate voltage is small enough that 
no electrons are essentially depleted. However, as men- 
tioned above, when the barrier hight is larger than the 
Fermi energy there exists no electrons to screen the ex- 
ternal potential and the potential distribution must be 
calculated sclf-consistently. 

Another adjustable parameter which can be accessed 
experimentally is the geometry of the structure. Of 
course, in the simplistic models describing QPCs this 
does not play an important role, since the constriction is 
assumed to be isotropic in the current direction, in con- 
trast to the experimental findings. It is well known that, 
the shape of the QPCs, as well as the cooling and biassing 
procedure, is important when measuring interference or 
magnetic focusing. In FigjSJ we compare two different 
gate patterns considering typical gate separations W for 
a fixed gate voltage, Vg = —2.0 V. The smooth configu- 
ration (CI), exhibits a stronger non- linear behavior when 
W is changed from 200 nm (black dashed line) to 300 nm 
(red dashed line). This relics on the fact that, at the first 
order approximation, the screening is better when more 
electrons are accumulated at the opening of the QPC. 
However, since the screened potential V SCI (x,y) can be 
obtained from the external potential V CKt (x,y) via the 
Thomas-Fermi screening, 

K q CI = VlJe(q), (1) 

where e(q) = 1 + 2/(a* B \q\), is the Thomas-Fermi dielec- 
tric function with q being the momentum and a B (= 9.8 
nm for GaAs) the effective Bohr radius, the long range 
fluctuations (large g), compared to a Bl are less screened 
whereas the short range fluctuations are predominantly 
screened. Considering the sharp configuration (C2) this 
observation becomes more evident, since the potential 
profile across the QPC varies smoother than of the con- 
figuration (CI), when varying W . From an experimental- 
ists point of view, therefore, drawing shaper QPC struc- 
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FIG. 3: (Color online) Potential profile across the QPCs for 
two different geometrical patterns. Insets depict the smooth 
(CI) and the sharp edged (C2) patterns. Dark (blue) regions 
are electron depleted, i.e. the local potential is larger than the 
Ef, denoted by the solid thick horizontal line. White contours 
denote the depletion boundary, where the width of the QPC 
is set to be 200 nm and the curvature is changed. The dashed 
lines stand for the first configuration and two different W 
values (= 200 nm (black) and 300 nm (red)), whereas for the 
sharper configuration four W values are selected. 



tures by electron beam lithography may lead to a better 
(linear) control of the potential profile which is closer to 
the ideal potential profile. This is certainly in contrast 
to what we would expect from an non-interacting model, 
however, it is known by the experimentalists that defin- 
ing the QPCs with sharper edges increase the quality of 
the visibility signaP^. For the C2 configuration we also 
observe that, the potential profile becomes almost insen- 
sitive to the width W above 300 nm, which coincides with 
our previous finding of better screening of the long range 
fluctuations. It is worth to note that, our calculation 
scheme is beyond the simple Thomas-Fermi screening 
scheme in obtaining the bare confinement potentiaP^. It 
fully takes into account the interaction effects, however, 
does not include any quantum mechanical effects. In a 
better approximation, of course, one should also solve the 
Schrodingcr equation self-consistcntly in 3D. This proce- 
dure is known to be costly in terms of computational 
cost even only if the 2DES is treated quantum mechani- 
cally. Since we are interested in either zero or very strong 
magnetic fields, representing electron as a point charge is 
still a reasonable and valid approximation. We will dis- 
cuss the justification of this assumption in the presence 
of an external magnetic field in more detail in Sec. |III| 
where we also discuss the limitations of our model. 

The different configurations of the gate patterns are 
also important when investigating the scattering pro- 
cesses by magnetic focusing experiments^. It is apparent 
that the scattering patterns of the electron waves will not 
only depend on the impurity distribution but also on the 



FIG. 4: Potential cross-section of (a) gated and respectively 
(b) etched QPCs at selected gate biases and etching depths 
for CI at a fixed W = 150 nm. Insets focus on the high bias 
or shallow etching profiles 



structure of the QPCs. We expect that for the sharper 
defined QPC (C2) the scattering should depend weakly 
on the QPC opening, since the width does not change 
along the constraint. Whereas, for CI configuration small 
changes at W should affect the scattered waves drasti- 
cally. Another comment on the experimental setups is to 
the RoddarcW experiments, since the formation of (sta- 
ble) integer filling region is important to explain^ the 
findings, we believe that C2 type configurations would 
be leading to a better resolution of the transmission am- 
plitudes. The choice of the structure and the width ap- 
parently depends on the experimental interest, which is 
believed to be irrelevant when modelling QPCs as ideal 
point contacts, and the QPCs are not only denned by 
gates but alternatively also by chemical etching. The 
gated structures are of course more controllable, how- 
ever, at high gate voltages required for depletion, electri- 
cal sparks can occur, therefore the structure can even be 
destroyed. In such situations etching denned QPCs are 
preferred, although without further gates one loses the 
full control of the potential profile. In the next section, 
we will compare the potential profiles of etched and gate 
defined QPCs, to show that in some cases etch defined 
QPCs may be more useful to obtain steeper potential 
profiles. 



B. Etching versus gating 

For simple calculation purposes, QPCs are modelled 
either as a finite potential well or with a parabolic con- 
finement potential, perpendicular to the current direc- 
tion. Starting from the early experiments 1 , usually the 
conductance is measured as a function of the applied gate 
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voltage which presents clear quantized values. This quan- 
tization can be well explained by the Landauer formula^ 
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where ballistic transport is assumed to take place, i.e. 
the transmission is given by |t n ,m| 2 = ^n,m, and no chan- 
nel mixing is allowed. The (integer) number of chan- 
nels N c is defined by the Fermi energy and the width 
of the constriction, in general. The gate defined QPCs, 
at a first order approximation, can be represented by 
parabolic or finite well potential profiles. However, it is 
known that the chemical etching process creates (side) 
surface charges, which in turn generates a steeper poten- 
tial profile at the edges of the sample. In this situation 
it is apparent that, the confinement potential can not 
be assumed to be parabolic, rather a steeper potential 
should be considered. In this section we compare these 
two different constriction profiles, namely the gated and 
the etched ones. 

The self-consistent potential across the QPC at the 
center is plotted versus the lateral coordinate in Fig. [4] 
The 2DES under the gates is depleted at the gate volt- 
ages larger than —1.5 V, similarly when the depth of the 
etching is larger than 47.7 nm. In both cases, till the de- 
pletion is reached the potential profile is varying rather 
smoothly and the depth of the potential depends linearly 
on the applied gate bias or the depth. This behavior 
is drastically changed when electrons are depleted, the 
potential now strongly depends on the gate potential or 
depth. Moreover, for the etched structure, potential be- 
comes very steep at the edge of the QPC when D c > 47.7 
nm, i.e. is deeper than the donor layer. The transition 
from linear behavior to non-linear behavior is simply due 
to the significant change of the screening properties of 
the system. Once the electrons are depleted, the exter- 
nal potential can no longer be screened, therefore the 
amplitude of the self-consistent potential increases by a 
large amount. Therefore, screening calculations based on 
the formula given in Eqn. can not account for such 
situations, where the dielectric function is not aware of 
the Fermi energy, i.e. the occupancy. A better approxi- 
mation to this approach is to consider a Linhard type di- 
electric function, which also takes into account the Fermi 
momentum^. It is known that, such an improvement 
will also cover some of the quantum mechanical aspects 
(such as the wave functions), which brings extra oscilla- 
tions to the potential profile 2 ^. However, for our present 
interest we neglect this correction knowing that the self- 
consistency of the calculation scheme already takes into 
account the occupation and the ID electron density at 
the QPC satisfy the validity condition n c \a* B >> 1 of the 

tfaeh. 

We summarize our findings in Fig. [5] where we show 
the electron density (left) and potential profiles (right) 
for typical gate biases V g and etching depths D c versus 
the spatial coordinate. We choose a representative cross- 



section of the obtained profiles along the current direc- 
tion x, where the y coordinate is fixed at 450 nm. Figure 
[5^ depicts the density profile for selected depths of etch- 
ing varying from shallow (D c — 4.7 — 19.1 nm) to deep 
(D c = 38.2 - 78.4 nm). For the smallest two D c 's, the 
2DES is not depleted beneath the pattern and the density 
profile is rather smooth. Depletion is observed when the 
depth is larger than 19.1 nm, however, until the etching 
depth reaches to the depth of 2DES (~ 60 nm) we do not 
see the surface charges (the spike like point, indicated by 
the arrow) at the level of the 2DES. The inset of Fig. 5jp 
shows the electron density distribution in a color coded 
contour plot together with the corresponding potential 
profile across the white (dashed) line. The thin (green) 
lines contouring the depleted (red) region indicates the 
spatial distribution of the surface charges. The potential 
is steeper compared to that of the gated one (Fig. [5ji) 
and the profile does not show any considerable variation 
once the etching depth reaches the plane of the 2DES. 
This behavior clearly exhibits the uncontrollability of the 
etched samples, since the corresponding potential profile 
obtained for the gated samples vary slowly on the length 
scales of the Fermi wavelength, even if the 2DES is com- 
pletely depleted beneath the gates. Moreover, the am- 
plitude of the potential strongly depends on the applied 
gate voltage. The slow variation of the potential is not 
the case for the etched sample, for example consider the 
case when D e is changed from 57.3 to 78.4 nm, and com- 
pare it with that of the gated sample when voltage is 
changed from —2.0 V to —3.0 V. For the gated sample 
the potential profile remains still smooth, however, the 
spatial distribution of the electron density is almost un- 
changed. Meanwhile, for the deeply etched sample both 
the electron density and the potential profile are steep 
and the steepness depends very weakly on the etching 
depth. 

From the above discussion we conclude that, for the 
gated samples the electron density distribution is weakly 
effected by the applied gate potential when the depletion 
is once obtained, meanwhile the potential is smooth and 
strongly depends on V s . For the etched samples, poten- 
tial profile becomes very steep when the etching depth 
exceeds the depth of the 2DES, since the (side) surface 
charges pin the Fermi level at the mesa surface to the 
mid gap of GaAs forming a Schottky like barrieil^. We 
should also note that, at zero bias, more electrons are 
populated under the gates, which is not the case for the 
etched samples. As a rule of a thumb, when a steeper 
potential profile is required one should consider chemi- 
cal etching where the etching depth is deeper than the 
electron layer and one should keep in mind that biasing 
gates with high potential does not necessarily imply that 
the electron density profile is also changed considerably. 

The outcome of the self-consistent solution of the 3D 
Poisson equation considering QPCs at zero magnetic field 
is two fold: i) the geometrical properties (i.e. consider- 
ing CI or C2 type patterns) strongly change the potential 
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FIG. 5: Spatial distribution of the electron density (a,c) and 
screened potential (b,d) for etched (upper panel) and gated 
samples (lower panel), y — 450 nm. 



landscape in the close vicinity of the QPCs. We found 
that the smoother constrictions with a larger width W 
can be modelled better with the ideal point contacts, i.e. 
parabolic confinement. The sharper constrictions can 
be considered as finite well profiles, up to a first order 
approximation and potential profile remains unchanged 
when considering W > 300 nm. ii) Due to surface pin- 
ning the etched samples, generate steeper potential pro- 
files and the density profile remains unaffected once etch- 
ing is deeper than the depth of the 2DES. These numer- 
ical results, show a strong deviation from the the widely 
used "ideal point contact" and "frozen charge" models, 
proposing that depending on the experimental interest it 
is important to reconsider the geometrical (CI, C2) and 
structural (gated/etched) factors defining the QPC un- 
der investigation. As final remark, the artifacts, such as 
local minima and maxima at the potential and density 
profiles, resulting from the previous non-sel f consist ent 
Thomas- FermP^T and "frozen charge" models^ 9 -* 31 * 36 * are 
resolved by considering the 3D calculation scheme. 



III. FINITE B FIELD 

The aim of this and the next section is to calculate, 
and compare, the density and subsequently the current 
distributions, in the close vicinity of the QPCs, within 
interacting and non-interacting models in the presence 
of an external perpendicular magnetic field. Here we 
take into account electron-electron interaction within 
the Thomas-Fermi theory of screening also considering 
strong magnetic fields using the potential profile calcu- 
lated in the previous section, as an initial configuration 
of the landscape. We compare the spatial distribution 
of the Landauer-Biittiker edge states (LB-ES) with the 
distribution of the incompressible edge states, where the 



applied external current is confined' 9 * 40 * 47 ' * 48 !. Next we 
discuss the limitations of the TEA, and suggest improve- 
ments on the calculation scheme based on i) quantum 
mechanical considerations, such as the finite extension 
of the wave functions, and ii) replacement of the global 
density of states with the local one. 

A. Thomas-Fermi-Approximation (TEA) 

The enormous variety of the theories describing the 
density and current di stributions at the quantum Hall 
systems^* 21 * 38 * 40 * 47 * 48 * 49 1 already show the challenge in giv- 
ing a proper prescription to these quantities. These the- 
ories can be grouped into two: the current is carried ei- 
ther (i) by the compr essible regio n d 38 * 50 * 5 - 1 ^ or (ii) by the 
incompressible regiond^* 40 * 47 * 48 * 5 ^ . Moreov er (an d confus- 
ingly) these regio ns can reside a t the bull j 11 * 49 ^ or at the 
edge of the sampl e 9 * 38 * 40 * 47 * 51 * 53 !, depending on the model 
considered and the magnetic field strengtlP. For the sake 
of completeness, we start with a generic Hamiltonian de- 
scribing an electron subject to high magnetic fields. 

H° = H Q + Ct + V ext + V£, (3) 

where cr(= ±1/2) is the spin degree of freedom, Ho the 
kinetic part, V cxt and Vi nt the external and the inter- 
action potentials, respectively, and V£ Zeeman terrrP^. 
Our first assumption is to neglect the spin degener- 
acy, knowing that the effective band g— factor for the 
GaAs/AlGaAs hetero-structures is a factor of four less 
than the one of a free electron gas, and therefore Zeeman 
splitting is much smaller than the Landau splitting (hw c , 
with uj c = eB/m*, i.e. \g* HBB\/hw c ~ 0.027, where hb 
is the Bohr magneton. However, Zeeman splitting can be 
as large as the Landau splitting if exchange and correla- 
tion effects are taken into account at significantly high 
magnetic fields, hence filling factor v (— Ep/huj c ) one 
plateau can be observed experimentally. On the other 
hand, for higher filling factors v > 1 exchange and corre- 
lation effects are assumed to be small and Zeeman split- 
ting is negligible. Thus one can consider spin-less elec- 
trons in the magnetic field interval we are interested in, 
which yields to the effective Hamiltonian, 

H cH = H + K x t(z, y) + VintO, v)- (4) 

The kinetic part, H , can be solved analytically using 
the Landau gauge which yield plane wave solutions in 
one direction (y) and Landau wavefunctions in the other 
direction (x). Here, implicitly, a long Hall bar (ideally in- 
finite) is assumed, which is justified while the Fermi wave 
length is (~ 30 — 40 nm) much smaller than the sample 
length under consideration (L y ~ 1500 nm). Then the 
eigenfunctions of Hq can be expressed as 

k (x, v) = — = exp (ik v .y) x 

' v y/2 n n\y/^lL y V 

exp[-(^) 2 /2]x J ff„(^), (5) 
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where k y is the quasi-continuous momentum in y direc- 
tion, n the Landau index, X = —l 2 k y a center coordi- 
nate, and H n (^) the nth order Hermite polynomial with 
the argument £, whereas the eigen energies are 

E n = hio c {n+ 1/2). (6) 

The essence of the TFA relies on the fact that the poten- 
tial profile varies smoothly on the quantum mechanical 
length scales. Through out this work we will only con- 
sider the 6 — 8 T interval and as a rough estimation the 
extend of the wavefunctions in x direction (of the ground 
state, i.e. v = 2), or the magnetic length / = ^h/muj c , 
will be always similar or less than 10 nm, therefore in 
almost all cases neglecting the finite extend of the wave- 
functions is still reasonable. However, we have already 
seen that for the etched samples the potential is quite 
steep and the results obtained from the TFA may be 
doubted, which we will address in the next section. 

At the moment let us consider a case that the con- 
dition of TFA holds, i.e. the total potential Vt t(x,y) = 
V ex t(x, y)+Vi n t(x, y) varies smoothly on the quantum me- 
chanical length scales and the sample is long enough (i.e. 
kpLy 3> 1). Then we can replace the wavefunctions in 
both directions by wave packets centered at X, and at y 
and the center coordinate dependent eigenenergy E n {X) 
can be approximated to E n + Vtot (X, y) . It follows that 
the spatial distribution of the el ectron d ensity within the 
TFA is given by the expression l 32 l 33 l 55 l ; 

n el (x,y) = J D{E,{x,y))f{E + V tot (x,y) - ^*)dE (7) 

with D(E,(x,y)) (local) density of states, f(E) — 
1 / '[exp(E I fc&T) + 1] the Fermi function, u* the electro- 
chemical potential, which is a constant in equilibrium 
state, fee the Boltzmann constant, and T the tempera- 
ture. Once the electron density is obtained, the interac- 
tion potential, i.e. the Hartree potential, can be obtained 
from 

2e 2 f 

Vint(x,y) = — / K(x,y,x',y')n(x',y')dx'dy'. (8) 

Here, R is an average dielectric constant (= 12.4 for 
GaAs) and K(x,y,x' ,y') is the solution of the 2D Pois- 
son equation satisfying the boundary conditions dictated 
by the sample. The results reported in this and the 
following section assume periodic boundary conditions, 
where a closed form of the kernel K(x,y,x' , ?/) can be 
obtained analyticallj^. Equations Q and (TsJ form a 
self-consistent loop to obtain the potential and the den- 
sity profiles of a 2DES subject to high perpendicular mag- 
netic field in the absence of an external current at equilib- 
rium, which has to be solved iteratively using numerical 
methods. The computational effort to calculate electron 
and potential profiles within the TFA is much less than 
that of the full quantum mechanical calculation proce- 
dures. The results of both agree quantitatively very well 
in certain magnetic field intervals where the widths of the 



incompressible strips Wis ( m which the potential changes 
strongly) is larger than I. If Wis ^ I condition is reached, 
first of all the TFA becomes invalid and the calculation 
of the electron density should include the finite extend of 
the wavefunctions. The underestimation of the quantum 
mechanical effects lead to existence of artificial incom- 
pressible strips both in the no n self- consistent electro- 
static appro ximati on (NSCESA/2ZES1 and self-consistent 
TFA schemes?™. In fact, as early as the NSCESA, the 
self-consistent schemes which also took into account finite 
extend of the wave functions already pointed out the sup- 
pression of the incompressible strips in certain magnetic 
field intervals^ and also in the recent reportsP^Ql. For a 
systematic comparison of the calculated widths of the in- 
compressible strips within the TFA and the full Hartree 
approximations, we suggest the reader to check Ref. [9 , 
where a simpler quasi-Hartree scheme is proposed to re- 
cover the artifacts arising from TFA. 



B. Corrections to the TFA 

Historically, the first implementation of the TFA, in- 
cluding electron interactions, to quantum Hall systems 
goes back to the seminal work by Chklovskii et.aP^. 
There it was shown that within the electrostatic approx- 
imation the 2DES is divided into two regions which have 
completely different screening properties. In this model, 
a translation invariance is assumed in the current (y— ) 
direction. Due to finite widths of the samples in the 
x direction the electrostatic potential is bent upwards 
at the edges of the sample, hence the Landau levels, 
E n (X) = E n + V(X). Inclusion of the Coulomb interac- 
tion and the pinning of the Fermi energy to the Landau 
levels result in two regions (strips): i) The Fermi energy 
is pinned to one of the highly degenerate Landau levels, 
then the screening is perfect, effective potential is com- 
pletely flat (metal like) and electron density varies spa- 
tially ii) The Fermi energy falls in between two consecu- 
tive Landau levels, screening is poor, effective (screened) 
potential varies (the amplitude of the variation is Hlo c ) 
and electron density is constant over this region. It is 
apparent that, if the potential varies rapidly on the scale 
of I, the TFA fails and the results become unreliable. 
This condition is realized when considering narrow in- 
compressible strips having a width smaller than the mag- 
netic length. Therefore, one should include the effect of 
wave functions within these narrow strips. One way is, 
of course, to do full Hartree calculations. We already 
mentioned the challenges in the computational effort. A 
simpler approach is to replace the delta wave functions of 
the TFA with the unperturbed Landau wave functions, 
i.e. quasi Hartree approximation 9 (QHA). The findings 
of the QHA is shown to be more reasonable than of the 
TFA, which now also includes the finite extend of the 
wave functions in the close vicinity of the incompressible 
strips. Therefore, as an end result, when the Wis ^5 I con- 
dition is reached the incompressible strip disappear due 
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to the overlap of the neighboring wave functions. Based 
on this fact, in our calculations in the following we will 
exclude the effects arising from the artifacts of the TFA 
by considering a spatial averaging of the electron den- 
sity on the length scales smaller than I, which is known 
to b e releva nt in simulating the quantum mechanical ef- 

fect d 41 | 5» | HU | 

We should also make one more point clear that with 
the NSCESA usually a gate defined quantum Hall bar 
is considered. It is more common to define Hall bars by 
chemical etching and the edge potential profile is much 
more steep compared to gated samples which was shown 
in the previous section. Therefore, to fit the predictions 
of this model concerning the widths and the positions of 
the incompressible strips with the experimental data one 
has to assume thatPQ i) the 2DES and the gates are on 
the same plane and ii) the gate voltage applied should 
be fixed to the half of the mid gap of the GaAs, i.e. pin- 
ning of the Fermi energy at the GaAs surface. In fact, 
after making these two crucial assumptions the experi- 
mental findings of E. Ahlswede^ perfectly fits with the 
NSCESA. However, the widths (and the existence) of the 
incompressible strips strongly deviate from the predic- 
tions, since only the innermost incompressible strip can 
be observed. We have argued that, the widths of the in- 
compressible strips strongly depend on the slope of the 
potential, i.e. if the external potential is steep the in- 
compressible strips are narrow. Therefore one can easily 
conclude that, since the widths of the outermost incom- 
pressible strips become smaller than the magnetic length, 
the outer incompressible strips, i.e. the ones close to the 
edge, disappear and could not be observed. The overesti- 
mation of the Wis within the TFA becomes more severe 
when an external current is imposed to the system, which 
we will discuss in the next section. Before discussing the 
results of the relaxed TFA, we want to touch another lo- 
cally defined quantity, namely the DOS, and comment on 
the implementation of the global DOS to our local TFA. 

In the absence of impurity scattering the DOS of an in- 
finite (spin-less) 2DES is given by the bare Landau DOS 
as 



D(E) 



1 

2ttP 



(9) 



however, this DOS is broadened by the scattering pro- 
cesses, whic h ca n be described in self-consistent Born ap- 
proximation^^ accurately for short range impurity po- 
tentials yielding a semi elliptic broadening. Of course, 
other impurity models and scattering processes can also 
be considered resulting in Gaussian or Lorentzian broad- 
ened Landau DOSSES, In such descriptions of the DOS 
broadening an infinite 2DES is assumed and the DOS 
is calculated for impurity distributions averaged over all 
possible configurations. Inserting this (global) DOS in 
Eqn. Q can be justified again if the TFA condition is 
satisfied. We have already shown that this condition is 
violated when a narrow incompressible strip is formed 
where the external potential is poorly screened. There- 



fore, the actual distribution of the impurity potential be- 
comes more effective at these transparent regions. We 
should also note that the effect of screening on the DOS 
of an infinite system has been investigated in detail in 
Ref. [S3] and it has been shown that, since the screen- 
ing is poor within the incompressible regions, the DOS 
becomes much broader than that of the non-interacting 
case. Hence, the gap between two consecutive Landau 
levels is narrower within the incompressible strips com- 
pared to compressible strips. Moreover, recently it has 
been shown that the (local) electric field within the sam- 
ple also leads to broadening of the (local) DOS^. The 
idea is basically that one calculates the Greens function 
for the given potential profile, which is a function of the 
applied magnetic field and external current, and obtains 
the local DOS from the general expression 

D(E, (x,y)) = J2 \$n,k y (x,y)\ 2 6(E-E n , ky ), (10) 



where <fr n ,k y (%, y) is the nth cigenfunction of the Hamil- 
tonian given at Eqn. (UJ with the eigenvalue E n .k y ■ In 
our above discussion about the formation of the com- 
pressible/incompressible strips we have mentioned that 
the potential varies locally whenever an incompressible 
strip is formed, where the variation is linear in position 
up to a first order approximation. Now let us consider 
a linear potential profile and re-obtain the local DOS 
(LDOS) following^ for the fcth Landau level, 



D k (E) = 



1 



2 fc+l fc l 7I -3/2|2 r 

with the level width parameter 

r = £j 



E H v \H k {E h /T)f (11) 



(12) 



where £ x — dV(x,y)/dx is the electric field in the x 
direction and E k = E - T 2 /(4fiLu c ) + (2k + l)huj c . The 
immediate consequence of a strong electric field in the x 
direction is a broadening of the LDOS, which happens 
at the incompressible strips. On the other hand, since 
£ x vanishes at the compressible strips, the bare Landau 
DOS is reconstructed from Eqn. ( 11 ) in the r — > limit. 



In summary, the TFA should be repaired when the 
widths of the incompressible strips is comparable or less 
than the magnetic field, since (i) the quantum mechan- 
ical wave functions have a finite extend and do overlap 
with the neighboring ones (ii) the LDOS are broadened 
where strong electric fields exist, (iii) within the poor 
screening regions the (Landau) gap is smaller than the 
ones of the nearly perfect screening regions. At the very 
simple approximation this artefact of TFA is cured by a 
spatial averaging over the Fermi wavelength resulting in 
nonexistence of narrow incompressible strips. In the fol- 
lowing an appropriate averaging process will be applied. 
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FIG. 6: The distribution of the spin degenerate LB-ES at 
(a) v « 4 and (b) v ~ 8. The potential cross-section at v = 8 
plateau (c). Sketch of the conductance (G) and expected co- 
herence S (d) for CI considering W = 150 nm with an applied 
gate potential V g = 2.0 V, so that the 2DES is depleted be- 
neath the gates. Since spin degeneracy is assumed each edge 
state carries two units of quantized current. 



C. Results 

The transport through the QPCs in the absence of an 
external magnetic field is well described by the Landauer 
formalism, summarized in Eqn. The main idea is 

that the transport is ballistic and due to the cancella- 
tion of the velocity and a ID D092I the conductance is 
integer multiples of the conductance unit, e 2 /h. These 
integers are just the number of channels N c , i.e. the num- 
ber of eigen energies below the Fermi energy. A similar 
path is followed when considering an external B field, 
which assumes that the conductance is ballistic within 
ID channels and neglects the electron interactions. These 
ID channels are formed whenever the Fermi energy co- 
incides with a Landau level (Landauer-Biittiker picture) . 
Before proceeding with the full self-consistent solution 
of the density and current distribution problem, we first 
investigate the positions of the Landauer-Biittiker edge 
states. The procedure is simple, to obtain the energy 
dispersion we use the relation 



E n ,k v (x, y) = E n + V{x, y) 



(13) 



and follow the equipotential (energy) lines coinciding 
with the Fermi energy. By doing so we would be able 
to discuss qualitatively the phase coherence S, taken to 
be equal to 1 when there is a full phase coherence and 
when there is none. We neglect all the external sources 
of decoherence and assume that the LB-ES are coherent 
at the length scales we are interested in, i.e. the wave 
functions of the associated channels do not overlap. 

In Fig. [6] we show the spatial distribution of the ex- 
pected positions of the LB-ES at two filling factors. The 



color scale depicts the self-consistent potential, whereas 
the black strips show the LB-ES. The white shaded ar- 
eas are the electron depleted regions. The v — 2 and 
v = 4 edge states nicely show the expected distribution 
which are spatially ~ 40 nm apart, Fig. [6^. Depending 
on the steepness of the potential or the magnetic field 
value, however, this distance may become less. For the 
etched samples (not shown) at the same filling factor the 
spatial distance between v = 2 and v = 4 edge states 
become almost half the value of the gated samples. For a 
filling factor v =8 plateau the outermost (the ones closest 
to the gates) two edge states are less than 15 nm apart 
from each other, and the wave functions extend over a 
larger distance. It is apparent that, when the two wave 
functions start to overlap, the coherence 5 is reduced 
drastically. However, for the ideal case (no overlap) H 
should stay constant for all plateau regions, since by def- 
inition the edge states can not cross each other. The 
conductance quantization is, of course, independent of 
the structure of the ES itself and according to Eqn. (|2| 
one should simply count the number of ES which cross 
the constriction. The conductance is shown by the sketch 
in Fig. [6]i, of course the sharp transition between the 
plateaus is changed when one considers level broadening 
or in general scattering. One should note that, although 
the LBES picture is useful in making qualitative argu- 
ments, one needs to grasp the actual distribution of the 
edge states to understand the physics observed at the 
experiments^. 

Next we investigate the distribution of the incom- 
pressible strips calculated self-consistently described by 
Eqns. (7j|8!. The conductance through the QPC can be 



rewritten 



G 



(14) 



as conjectured in Ref. |27j . where f CO ntcr is the filling fac- 
tor at the very center of the QPC. It is apparent that, if 
this value is an integer, i.e. incompressible, the conduc- 
tance is quantized. Therefore, it is important to study 
this condition for a realistic QPC geometry. In the fol- 
lowing we first calculate the filling factor distribution in 
the absence of an external current and then obtain the 
current distribution in the next section. 

In order to cure the artifacts arising from TFA i) we 
consider a DOS broadened by a GaussiarP%iven by 



D(E) 



1 



^ eM-[En-E]yT? mp ) 
27T p 7^ (15) 



n=0 



imp 



with the impurity parameter r; mp , which is chosen large 
enough Ti mp /huj c = 0.3 to cover the level broadening 
generated by the local electric field^ and self-consistent 
broadening effects^ and ii) a spatial averaging is car- 
ried out over the Fermi wavelength (~ 30 nm). Fig. [7] 
summarizes our results showing the spatial distribution 
of the incompressible ES, considering the quantum Hall 
plateau v = 2. Pedagogically, starting our investigation 
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from large magnetic fields is preferable; at large magnetic 
fields (Fig. [7k ), the system is mostly compressible (col- 
ored area) and the two incompressible (white) regions do 
not merge at the opening of the QPC. Therefore, both the 
Hall resistance and the conductance through the QPC is 
not quantized. As soon as one enters to the QH plateau 
almost all of the sample becomes incompressible shown 
in Fig. [7jl (in the absence of short range impurities) and 
both Rh and G becomes quantized. Decreasing the mag- 
netic field creates two incompressible ES which are spa- 
tially separated seen in Fig. [TJd, however, the quantiza- 
tion is not affected. At a lower magnetic field value these 
IS-ES disappear (Fig. [7k) as an end result of level broad- 
ening and (simulation) of the finite extend of the wave 
functions, now we are out of the QH plateau and G is 
no longer quantized. This picture and the LB-ES picture 
yield same behavior for the Rh and G, however, in the 
later one current is carried by the IS-ES, which we will 
discuss in the next section. The qualitative difference 
between the two pictures is the coherence as shown by 
the (red) dashed line in Fig. [7fc. First S presents min- 
ima in between two plateau regimes, since the IS die out, 
second at the higher edge of the QH plateau the system 
becomes completely incompressible therefore, it is not 
possible to define separate ESs hence coherence is lost 
(averaged). We believe that, this nonuniform behavior 
of the coherence within the QH plateau coincides with 
the experimental findings of Roche et.alpu however, we 
admit that other explanations are also possible. Another 
interesting experimental work is carried by the Regens- 
burg group, where they have investigated the amplitude 
of the visibility oscillations as a function of B field 18 . 
They have reported a maximum visibility at v « 1.5, 
which seems q uite o pposite to, what has been reported 
by other groups^EH. However, it is easy to see that their 
sample has a homogeneous width all over, which is not 
the case for other groups. From self-consistent calcula- 
tions^ it known that, if the sample width is narrower 
than 5-6 /im the center electron density (or filling factor) 
can differ strongly from that of the average one(s). An 
indication of such a case is also shown by numerical sim- 
ulations^. In interconnecting the Hall plateaus and the 
spatial distribution of the incompressible strips, we have 
used the findings of Ref. [9] where the current is shown 
to be flowing through the incompressible strips. This is 
in contrast to some of the modelsSSMl m th e literature 
where the opposite is proposed. In the next section, we 
will present the general concepts of the local Ohm's law 
and based on the absence of back scattering in the in- 
compressible strips we will show that the local resistivity 
vanishes and the external current should be confined to 
these regions. 
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FIG. 7: Spatial distribution if the incompressible strips 
(white areas) for characteristic B values (a) 6.8 T, (b) 7.3 T, 
(d) 8.3 T (e) 8.8 T calculated at temperatures hto c /k B T < 1, 
together with the sketch of Hall resistance and the coherence 
H (c). The QPC configuration considered here is CI with 
W — 150 nm, the gate voltage is chosen such that all the 
electrons beneath the gates are depleted. 



IV. 



CURRENT DISTRIBUTION WITHIN 
LOCAL OHM'S LAW 



The local (potential) probe experiment d 61 * 67 *^, 
brought novel information concerning the Hall potential 
distribution over the sample. The first set of experiments 
show clearly that the potential, therefore the current, 
distribution is strongly magnetic field dependent. It 
was shown that, out of the QH plateau regime the Hall 
potential varies linearly (Type I) all over the sample, a 
similar behavior to classical (Drude) result. Whereas, at 
the lower edge of the QH plateau the current is confined 
to the edges of the sample (Type II), which was shown 
to be coinciding with the positions of the incompressible 
strips. The most interesting case is observed when 
an exact (even) integer filling is approached. In these 
magnetic field values, the potential exhibits a strong 
nonlinear variation all over the sample, which was at- 
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tributed to the existence of a large (bulk) incompressible 
region. The explanation of these measurements acquired 
a local transport theory, where the conductivities and 
therefore current distribution can be provided also 
taking into acc ount in teractions. In the subsequent 
theoretical work d 9 ' 40 ' 69 ' the required conditions were sat- 
isfied and an excellent agreement with the experiments 
were obtained^. In the second set of experiments^ 
a single electron transistor has been placed on top of 
the 2DES and the local transparency, i.e. whether 
the system is compressible or incompressible, and the 
local resistivity have been measured. Comparing the 
transparency and the longitudinal resistivity, it has been 
concluded that the resistivity vanishes when the system 
is incompressible. 

Theoretically, if the local electrostatic potential and 
the resistivity tensor p(x,y) are known the current dis- 
tribution j(x,y) can be obtained from the local version 
of the Ohm's law 



E{x,y) = p(x,y).j(x,y) 



provided that 



E(x,y) = Vp,*(x,y)/e 



(16) 



(17) 



where the electrochemical potential is position dependent 
when an fixed external current / = J A j(x,y)dxdy is im- 
posed. In our calculations we assume a local equilibrium 
in order to describe the stationary non-equilibrium state 
generated by the imposed current, starting from a ther- 
mal equilibrium state obtained from the modified TFA. 
At this point if the local resistivity tensor is known, Eqn.s 
(16l-(17l should be solved once again iteratively for a 



given electron density and potential profile, where the 
equation of continuity 
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FIG. 8: The local current density calculated at different field 
strengths, same as in Fig[7] The intensity of the current den- 
sity is chosen such that, the applied current does not effect 
the density distribution, i.e. |j(a;,J/)| ~ 0.4 x 10 -3 A/m. 



V-j(r) =0 



(18) 



also holds. We assume that the local resistivity is related 
to the local electron density via the conductivity tensor, 
i.e. p(x,y) = (T _1 (nei(x, y)). For a Gaussian broadened 
DOS the longitudinal component of the conductivity ten- 
sor is obtained from 



2e^ 
h 



dE[- 



dE i 



OO 

£ 

n=0 



(n 



whereas Hall component is simply 



2e 2 



a H = 



(19) 



(20) 



where we ignored the self-energy corrections depending 
on the type of the impurity scatterers. We should em- 
phasize that , the above conductivities are used for consis- 
tency reasons, in principle, any other reaso nable impurity 
model like the commonly used SCBA'^EIl can be consid- 
ered. Assuming that TFA is valid, we can replace the 
local conductivities with the above defined global ones. 



In the absence of an external current our calculation 
scheme is as follows, we initialize Eqn. ^ using the to- 
tal potential obtained from the 3D calculations and ob- 
tain the electron density at relatively high temperatures 
(kT/huj c ~ 0.5) and use this density distribution to ob- 
tain resulting potential from Eqn. Next we keep on 
iterating until a numerical convergence is reached, where 
the electron density is kept constant. This is followed 
by the step where the temperature is lowered by a small 
amount and iteration process is repeated till the target 
temperature is reached. After the thermal equilibrium is 
obtained, we impose a small amount of external current 
and solve Eqns. ( 16|[i~7 1 self-consistently. While doing 
this second iteration, we fix the constant arising from the 
integral equation to a value such that the total number 
of electrons is kept constant with and without current. 
As a numerical remark, if the current loop does not con- 
verge we increase the temperature by a relevant amount 
and start the iteration procedure. The whole calculation 
scheme is composed of three different codes, which are 
written in C++, Fortran and Matlab respectively. In 
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order to obtain reasonable grid resolution parallel com- 
putation techniques were usecP^. 

Since we are interested in current distribution and also 
its effect on the density distribution we find appropriate 
to present our results in two separate sections i) where 
the applied current is weak enough that the electron and 
potential distribution is unaffected, linear response, ii) 
the applied current is sufficiently large so that the im- 
posed current induces a considerable change on the po- 
sition dependent electrochemical potential, out of linear 
response. 



A. Linear response regime 

The crucial part of the local approach is that for a given 
(large) magnetic field we can calculate the local poten- 
tials and electron distributions self-consistently. The re- 
sult of such a calculation is that the 2DES is essentially 
separated in two regions, i.e. compressible and incom- 
pressible, therefore for a given (obtained) density we can 
calculate the local conductivities via Eqns. (ppOjl. Let 
us now discuss the distinguishing conductance properties 
of these two regions starting from a compressible region. 
At a compressible region the Fermi energy is pinned to 
one of the Landau levels, screening is nearly perfect, self- 
consistent potential is flat and filling factor is locally a 
non-integer. According to Eqn.(20l, the Hall conductiv- 



ity is a non-integer and the longitudinal conductivity is 
non-zero meaning finite back scattering. Now the clas- 
sically defined drift velocity, and also its quantum me- 
chanical counter part, is proportional to the electric field 
perpendicular to the current direction. We have seen that 
at the compressible region the potential perpendicular to 
the current direction is flat, therefore the x component of 
the electric field is zero, hence the drift velocity. Mean- 
while, at an incompressible region the Fermi energy is in 
between two Landau levels, the filling factor is fixed to 
an integer value and potential presents a variation per- 
pendicular to the current direction. Due to the Landau 
gap the longitudinal conductivity vanishes, whereas the 
Hall conductivity assumes its (quantized) integer value. 
If one calculates the inverse of the conductivity tensor for 
the longitudinal component one obtains, 



Pi{x,y) 



*(x,y) + a 2 H (x,y) 



(21) 



thus the longitudinal resistivity vanishes within the in- 
compressible region pointing the absence of back scatter- 
ing. Of course, the simultaneous vanishing of both the 
longitudinal resistivity and the conductivity is a result 
of applied external (and perpendicular) B field and is 
obtained only in two dimensions. Moreover, since the x 
component of the electric field is now non-zero, the drift 
velocity is finite and the current is confined to this region. 
Combining these two one concludes that, if there exists 
an incompressible region somewhere in the sample all the 
external current is confined to this region otherwise (if 



there are no incompressible regions and all the system 
is compressible) the current is distributed according to 
Drude formalism, i.e. the current density is proportional 
to the electron density. 

We start our discussion of the current distribution 
when a small current is imposed for which the electro- 
static and electrochemical potential satisfies the linear 
response relation 

V((x, y);I) - V((x, y); 0) « »*((x, y);I) - (i* eq . (22) 

This condition essentially states that, the imposed cur- 
rent does not modify the electrochemical potential there- 
fore the electron density remains unchanged. Fig. [8] 
presents the current distribution which is calculated for 
the density distribution shown in Fig. [7] The correspon- 
dence between the positions of the incompressible strips 
and the current density is one to one. In the out of 
plateau regime the current is essentially distributed all 
over the sample, where no incompressible regions exist, 
Fig|8ji. As soon as one enters the QHP, i.e. when a large 
bulk incompressible strip (region) is formed, the essen- 
tial future of current distribution is not effected strongly, 
however, in this situation current is flowing in the in- 
compressible region. Tracking the positions of the ISs in 
Fig. [7J), we can readily guess the distribution of the cur- 
rent density in Fig. [SJj. Following our arguments about 
the smearing out of the narrow incompressible strips, we 
have a situation in which, again, the current is spread 
over the sample shown in Fig. [8ja,. Although, the IS 
vanishes the reminiscence of it still provides a narrow 
strip of small longitudinal resistivity and, therefore a 
higher amount of current is kept confined to these re- 
gions. Fig. [§];, presents the corresponding longitudinal 
resistance, when measured as a function of B together 
with the conductance across the QPC. The relation be- 
tween G and Rl is interesting, the conductance is quan- 
tized, as soon as Rl vanishes at large fields, however, 
becomes non-quantized even though the Rl = at the 
lower edge of the zero resistance state. Let us first discuss 
the B dependence of the Rl , it is finite if the system is 
compressible and is zero if an incompressible strip perco- 
lates from one edge of the sample to the other edge in the 
current direction, i.e. from source to drain. Therefore, 
existence of an IS percolating is sufficient enough to mea- 
sure zero longitudinal resistance. However, to have a con- 
ductance quantization the center of the QPC should be 
incompressible, which is a stronger restrictiorP-H Hence, 
in the lower edge of the QHP, the IS percolates but the 
center of the QPC is compressible. The implication of 
this fact to the coherence is a bit more complicated, we 
have seen that as soon as one enters to the QHP a large 
bulk IS is formed therefore the phase of the electrons 
is highly averaged. This implies that the coherence is 
relatively less than that of the two well separated ISs. 
On the other hand, at the lower edge of the conductance 
plateau, the ISs become narrower and are less immune to 
decoherence effects arising from the environment, hence, 
the coherence is reduced. Our above discussion coincides 
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with the recent experiments performed in small Mach- 
Zchnder interference devices (MZI), where the visibility 
is measured as a function of the external magnetic fields 
It is fair to note that, some other mechanisms providing 
B dependence of decoherence can also account for such 
a behavior. 

In the mentioned MZI experiment d 12 * 15 * and also the 
measurements performed at the group of S. RoddarcS^ 
a finite (and large) source drain voltage Vsd is applied 
either to measure the Vsd dependency of the visibility or 
the transmission. The intensity of the applied current, 
in these experiments, can not be treated within the lin- 
ear response regime, where the electrochemical potential 
remains constant, i.e. position independent. In the next 
section, we present the current and the density distribu- 
tion calculated where Eqn. (22 1 does not hold any more. 



B. Beyond linear response 

In the absence of an external current an equilibrium 
state is obtained by solving the Eqns. ^ and ^ self- 
consistcntly. Even in the presence of a small current, 
a Hall potential is generated which, in principle, modi- 
fies the electrochemical potential, i.e. tilts the Landau 
levels. This modification can be compensated by the re- 
distribution of the electrons, which certainly modifies the 
total electrostatic potential. If the applied current is suf- 
ficiently small, the modification is negligible, i.e. linear 
response. However, if the current is large, the resulting 
Hall potential is also large and one should re-calculate 
the electron density, and therefore the potential distri- 
bution till a steady state is reached. In this section we 
present the current and density distribution in the pres- 
ence of a large external current, where a local equilibrium 
is assumed implicitly. Fig. [9^ shows the electron density 
distribution in color scale for B = 7.1 T, where the inten- 
sity of the applied current is in the out of linear response 
regime. Note that the B value is chosen such that the 
Hall resistance is quantized, however, the conductance 
is not. The general behavior is similar to that of lin- 
ear response, however, it is clearly seen that the widths 
of the ISs are asymmetric with respect to y = 750 nm 
line, where current is driven in y direction. The asym- 
metry is induced by the large current, since the electrons 
are redistributed according to the new (self-consistent) 
potential distribution. The corresponding current den- 
sity distribution is plotted in Fig. |9]d, once more the one 
to one correspondence between the positions of the ISs 
and the local current maxima is apparent. The conse- 
quence of the asymmetry and thereby the widening of 
the ISs can be observed in the conductance and the Rl, 
such that the narrow IS at the upper side is smeared out 
much earlier than the one on the lower side. We should 
note that, such large currents heat the sample therefore 
the local temperature within the ISs is larger compared 
to the lattice temperature due to Joule heating^!. Such 
a (local) temperature dependence is treated explicitly by 
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FIG. 9: The local filling factor distribution and the corre- 
sponding local current density calculated at the default tem- 
perature. The current intensity is sufficiently high (\j(x, y)\ = 
2.0 x 10 -2 A/m) to induce an asymmetry on the density dis- 
tribution via position dependent electrochemical potential. 



H. AkerJ^l and his co-workers and a strong evidence is 
provided towards explaining the breakdown of the IQHE 
within this approach. Our present approach lacks such a 
treatment, therefore the competition between the widen- 
ing of the ISs and heating due to large currents is thought 
to be more complicated than presented here. As a con- 
sequence, the discussion of the coherence is far beyond 
our model, however, we think that the amplitude of the 
current when measuring visibility^ is assumably still in 
the linear response regime. 

In conclusion, by exploiting the local equilibrium and 
the properties of a steady state we have calculated the 
current distribution near a QPC in the out of linear re- 
sponse regime. We have shown that a asymmetry is in- 
duced on the density profile due to the bending of the 
Landau levels generated by the large Hall potential. We 
estimate that, the system can still be considered in the 
linear response regime if the ID current density is smaller 
than 0.042 x 10~ 2 A/m, which certainly depends on the 
details of the sample geometry. 
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V. SUMMARY 

In this paper, we provided a self-consistent scheme to 
obtain the electron density, potential profile and current 
distribution in the close vicinity of a QPC, within the 
Thomas-Fermi approximation. Starting from a litho- 
graphically defined 3D sample, we calculated the charge 
distribution at the surface gates, at the plane of 2DES 
and for etched samples at the side surfaces. The 3D self- 
consistent solution of the Poisson equation enabled us to 
present the similarities and differences between an etched 
and gate defined QPC. We found that, the relatively deep 
etched samples present a sharp potential profile near the 
edges of the sample. If the depth of the etching exceeds 
the depth of the 2DES, surface charges are calculated 
explicitly In the presence of a quantizing perpendicular 
magnetic field, we have calculated the distribution of the 
incompressible strips as a function of the field strength. 
We have argued that, if an incompressible strip becomes 
narrower than the magnetic length and/or if the trans- 
verse electric field is sufficiently large, due to Level broad- 
ening, the narrow incompressible strip is smeared. In the 
next step, the current distribution is obtained both in 
the linear response and out of linear response regimes 
using a local version of the Ohm's law assuming a steady 
state at local equilibrium. It is shown that the current is 
confined to the incompressible strips, due to the absence 
of back scattering, otherwise is distributed all over the 
sample. We have commented on the relation between 
the existence and percolation properties of the incom- 
pressible strips and the measured quantities such as the 
longitudinal resistance and conductance across the QPC. 
For the ideal clean sample, i.e. in the absence of long 
range fluctuations, it is shown that the QH plateau ex- 
tends wider than that of the conductance plateau. In the 
out of linear response regime, a current induced density 
asymmetry is presented for the first time in such geome- 



tries under quantized Hall conditions. The observable 
effects of such an asymmetry are not clarified, since it is 
also know that large currents increase the temperature 
locally due to Joule heating. 

A natural extension of the existing model is to include 
the spin degree of freedom and thereby exchange and cor- 
relation effect s^. A local spin density functional theory 
approach! 36 ! 37 ! is the much promising one among others 
such as Monte CarlrJ^ and exact diagonalization, from 
computational and application point of view. In fact, 
such an approach already exists, however, the current is 
handled within the Landauer-Biittiker formalism, which 
we think is not reasonable in the presence of large incom- 
pressible strips. On the other hand, a time dependent 
spin density functional model would be a good candi- 
date to describe current for the geometries under inves- 
tigation. The implementation of the Akera's theory, i.e. 
Joule heating, to our model is of course desirable which 
is been worked presently. 

Finally, in order to have a predictive power on the in- 
terference experiments we would like to utilize the exist- 
ing coherent transport models in describing the current 
together with our electrostatic model, which we are not 
able at the present. Another challenge is to simulate 
the real experimental geometries which already includes 
more then a single QPC and contacts etc. The numerical 
routine we developed is now able to do such large scale 
calculations within the linear response regime, however, 
still lacks describing the exchange and correlation effects. 
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